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Abstract 

The flow of viscous incompressible fluid over a periodically corrugated surface is investigated nu- 
merically by solving the Navier-Stokes equation with the local slip and no-slip boundary conditions. 
We consider the effective slip length which is defined with respect to the level of the mean height of 
the surface roughness. With increasing corrugation amplitude the effective no-slip boundary plane 
is shifted towards the bulk of the fluid, which implies a negative effective slip length. The analysis 
of the wall shear stress indicates that a flow circulation is developed in the grooves of the rough 
surface provided that the local boundary condition is no-slip. By applying a local slip boundary 
condition, the center of the vortex is displaced towards the bottom the grooves and the effective 
slip length increases. When the intrinsic slip length is larger than the corrugation amplitude, the 
flow streamlines near the surface are deformed to follow the boundary curvature, the vortex van- 
ishes, and the effective slip length saturates to a constant value. Inertial effects promote vortex 
flow formation in the grooves and reduce the effective slip length. 

PACS numbers: 68.08.-p, 83.50.Rp, 47.61.-k, 02.70.Dh, 47.15.Rq, 47.32.Ff 
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I. INTRODUCTION 



An accurate flow prediction in micro-channels is important for the optimal design and 
fabrication of microfluidic devices whose applications range from medicine to biotechnol- 
ogy na. 

The boundary conditions and the surface topology are major factors affecting 
the flow pattern near the solid boundary and the formation of recirculation zones. The flow 
separation at rough surfaces can modify the wall shear stress distribution or initiate insta- 
bility towards turbulence. In microfluidic channels, the vortex flow enhances the mixing 
efficiency 

a a 

and promotes convective heat transfer . In cardiovascular systems, 

the separation region at the entrance of branching vessels may trap lipid particles which 



could lead to arterial diseases 



10] . In the present study we examine the role of slip 



boundary condition in determining the flow properties near rough surfaces including the 
separation phenomena, and distribution of pressure and shear stress along the surface. 
Although the validity of the no-slip boundary condition is well accepted at the macro- 



scopic level^jrecent exp eriments 11 



lations 



HHQ, 



12 



20 



21 



3, 



14, 



15 



16] and molecular dynamics (MD) simu- 



221 ] reported the existence of a boundary slip in microflows. The 



model first proposed by Navier relates the slip velocity to the rate of shear via the propor- 
tionality coefficient, the so-called slip length. The MD simulations are particularly suitable 
for examining the influence of molecular parameters on the microscopic slip length at the 
liquid/solid interface. The advantage of the MD simulations is that a detailed flow analysis 
can be performed at the molecular level while the explicit specification of the boundary 
conditions is not required. In contrast to description of the flow near boundary by means of 
microscopic slip length, it is convenient to characterize the flow over macroscopically rough 
surfaces by the effective slip length, which is defined as a distance from the level of the mean 
height of the surface roughness to the point where linearly extrapolated bulk velocity profile 
vanishes. Recent MD studies have demonstrated that the effective slip length in a flow of 
and polymer melts 24j over a wavy surface agrees well with hydrodynamic 



simple fluids 



predictions [25 



23 



261 ] when the corrugation wavelength is larger than approximately thirty 



molecular diameters. 

The influence of surface roughness on fluid flow with either local no-slip or zero shear stress 



i.e. perfect s 



28|, 



29, 



30 



31 



ip) boundary conditions has been extensively studied in the past decades 



32 



27|, 



331 ] . see also a review section in 23]]. Analytical calculations have shown 
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that in a shear flow over a corrugated surface with microscopic no-slip or zero shear stress 



conditions, the effective boundary slip is insignificant macroscopically 



27 



29|. The effective 



no-slip boundary plane is located at the intermediate position between crests and valleys of 



the rough surface when the no-slip condition is imposed along the solid boundary 



31 



34 



16 



35| . For an arbitrary surface roughness with small amplitudes, the slip coefficient 



in the Navier model is proportional to the average amplitude of the wall roughness and 
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Applying the no-slip 



depends on the position of the origin of the coordinate system 
boundary condition along the wavy surface, Tuck and Kouzoubov [3jJ have demonstrated 
that the effective slip length is inversely proportional to the corrugation wavelength and 
quadratically proportional to the amplitude of the surface roughness. However, the series 



expansion method used in 



3JJ fails at large wavenumbers, ka > 0.5, when a backflow appears 



2J, |28|, [33j. Despite 



inside the grooves of the substrate. The effective slip length for a flow above the _surface with 
deep corrugations only weakly depends on the depth of the grooves 
considerable analytical efforts, the relation between the vortex flow structure in deep grooves 
and the effective slip length has not yet been systematically investigated. 

The laminar flow separation at the corrugated surface with the local no-slip bound ary 



conditions depends on the depth of the grooves and the Reynolds number 



40 



41 



36. 



37J, 



38 



39, 



42j, |43j. In a creeping flow over a sinusoidal surface, the flow circulation appears 



in sufficiently deep grooves and, as the corrugation amplitude increases, the vortex grows 



and remains symmetric 4jJ, |4J, |45(. With increasing Reynolds number, the vortex flow 



forms even in shallow grooves, the circulation region expands, and the center of vorticity 



is displaced upstream 36j, |39_|, |40|, |42| . In the limit of small-scale surface roughness and for 
no-slip boundary conditions, the apparent slip velocity at the mean surface becomes more 
negative as the Reynolds number increases [3l|. A noticeable change in the effective slip 
length was also observed at Re > 100 for laminar flow over deep grooves when the local slip 
length is comparable to the corrugation amplitude 24J. However, the influence of the local 
slip condition at the curved boundary on the vortex flow formation has not been considered 
at finite Reynolds numbers. 

This paper is focused on investigation of the effects of local slip boundary conditions and 
the Reynolds number on the flow structure near periodically corrugated surfaces and the 
effective slip length. We will show that for the Stokes flow with the local no-slip bound- 
ary condition, the effective slip length decreases with increasing corrugation amplitude and 
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the flow circulation develops in sufficiently deep grooves. In the presence of the local slip 
boundary condition along the rough surface, the effective slip length increases and the size of 
the vortex is reduced but its structure remains symmetrical. The analysis of the numerical 
solution of the Navier-Stokes equation with the local slip condition shows that the inertial 
effects promote the asymmetric vortex flow formation and reduce the effective slip length. 

This paper is organized as follows. The details of a continuum model and the imple- 
mentation of the local slip boundary conditions are described in the next section. The 
analytical results for the Stokes flow over a wavy surface by Panzer et al. 26] are summa- 
rized in Sec. lIII Al The analysis of the effective slip length and the flow structure is presented 
in Sec. lIII Bl for the no-slip case and in Sec. lIII CI for the finite microscopic slip. The effect of 
Reynolds number on the effective slip flow over a periodically corrugated surface is studied 
in Sec. lIII Dl A brief summary is given in the last section. 



II. DETAILS OF NUMERICAL SIMULATIONS 

The two-dimensional incompressible and steady Navier-Stokes (NS) equation is solved 
using the finite element method. The computational setup consists of a viscous fluid con- 
fined between an upper flat wall and a lower sinusoidal wall (see Fig. [I]). The corrugation 
wavelength of the lower wall is set to A and is equal to the system size in the x direction. 
The upper wall is located at h — A above the reference line at a = 0, which is defined as the 
mean height of the surface roughness. The continuity and NS equations are 

V-u = 0, (1) 

p (u ■ Vu) = -Vp + /iV 2 u, (2) 

where u = u i + v j is the velocity vector in the cartesian coordinate system, p is the pressure 
field, p is the fluid density, and /i is the Newtonian viscosity. 

The penalty formulation is employed to avoid decoupling between the pressure and ve- 
locity fields [461 ] . In this method, the continuity equation is replaced with a perturbed 
equation 

V-u = -f, (3) 
where A is the penalty parameter, which ensures the incompressibility condition. Thus, the 
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modified momentum equations in the x and z directions are 



p (u ■ Vu) = AV(V ■ u) + pV 2 u, 



(4) 



p (u • Vu) = AV(V • u) + /iV 2 f . (5) 

The advantage of the penalty formulation is the elimination of pressure and the continuity 
equation. The penalty parameter A must be large enough so that compressibility errors are 
minimal. The upper bound of A is determined from the condition that the viscous effects 

nn 

are not destroyed by the machine precision |46|, |47J. The penalty parameter A should be 
chosen according to the rule 

A = cmax{/i, pRe}, (6) 

where Re is the Reynolds number and the constant c is recommended to be about 10 7 for 
computations with double-precision 64 bit words [471 ] . 

The Galerkin formulation of Eq. (j3j) and Eq. ffih can be explicitly rewritten as 



Jn 



3N 3 



+ ViUj- 



dN 



Oz 



+ 



dNi fdNj dN , 
ox V ox oz 



+ 



Jn 



V I "a a + "IS » 



(7) 



P ^ UiVj ~d^ + ViVj ~dT 



+ 



Oz \ ox oz ' 



+ 



/dNi dNj dNdN,^ , 
V ax ax az a^ 



HffS'*, (8) 



where Afj is the weight function, Nj, Uj, Vj are the node shape function and the velocities 
in each element, U{ and Vi are the lagged velocities, and the right hand side (RHS X , RHS Z ) 
terms include the boundary velocities. 

In our simulation, the bilinear quadrilateral elements (i, j = 1, 2, 3, 4) with non-orthogonal 
edges are transformed to straight-sided orthogonal elements by introducing the natural co- 
ordinates £ = £(x,z) and rj = r](x,z). The shape functions N in the natural coordinate 
system are defined as 

Ni = (i+MKilM) j = 1 _ A (9) 
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where £j and ^ are the corner points of each element (see Fig. [2]). 
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In the next step, Eq. ((Tj) and Eq. (jSj) are integrated numerically using four-point Gaussian 



quadrature 
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The final system of equations is constructed as follows: 





" u ' 




RH S x 




. V - 




. RH S z . 



(10) 



where the terms RHS X and RHS Z contain the velocities at the boundary nodes. 

The boundary conditions must be specified at the inlet, outlet and upper and lower walls 
of the Couette cell. The periodic boundary conditions are imposed at inlet and outlet along 
the x direction. A finite slip is allowed along the lower wall while the boundary condition at 
the upper wall is always no-slip. In the local coordinate system (spanned by the tangential 
t and normal n vectors), the fluid velocity along the lower wavy wall is computed from 



u t = L Q [(n ■ V)«t + u t /R{x)), 



(11) 



where u t is the tangential component of u = u t t + u n n, L is the intrinsic (or microscopic) 
slip length at the flat surface, and R(x) is the local radius of curvature [26| . The radius of 
curvature is positive for concave and negative for convex regions. The Navier slip condition 
for a flat wall is recovered from Eq. (TTTT) when R(x) — > oo. The effective slip length L e g at 
the corrugated lower wall is obtained by extrapolating the linear part of the velocity profile 
(0.45 ^ z/h ^ 0.9) to zero velocity with respect to the reference line a = 0. 

The simulation begins by setting the no-slip boundary condition at the upper and lower 
walls as an initial guess. Once Eq. (TTUj) is solved, the fluid velocities at the lower boundary 
are updated using Eq. (fTTl) . This iteration is repeated until the solution is converged to a 
desired accuracy. The convergence rate of the solution remains under control by using the 
under-relaxation value 0.001 for the boundary nodes. The results presented in this paper are 
obtained with the grid resolution 150 x 150 in the x and z directions, respectively. In order 
to check the accuracy of the results, several sets of simulations were also carried out with 
a finer grid 180 x 180. The maximum relative error of the effective slip length due to the 
grid size is L e g/h = 0.003. The converged solution of the Navier-Stokes equation satisfies 
the following boundary condition: 

du t 1 1 1 

= T 



Ut 



LJx) 



dn 



LJx) 



R(x) 



(12) 



where L s (x) is the local slip length in the presence of surface curvature 26|. 
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The accuracy of the numerical solution is checked by the normalized average error, which 
is defined as 



[e i K (13) 



error 

' ^~~f I u 'i 

where iVp is the total number of computational nodes, w" is the velocity at the node i and 
time step n, and is the velocity in the next time step. The typical value of the error 
in the converged solution is less than 10~ 9 . Throughout the study, the results are presented 
in the non-dimensional form. The length scale, shear rate, shear stress, and velocity are 
normalized by h, 7*, t*, and U*, respectively, where 7* is the shear rate in the case of 
no-slip boundary condition at the flat upper and lower walls, and = /ry* and U* = hj*. 

III. RESULTS 

A. Analytical solution of the Stokes equation for viscous flow over a wavy wall 



The effect of small periodic surface roughness on the effective slip lengt 



1 has been previ- 



ously investigated for pressure-driven flows in a Couette geometry [25j, |26[ . The analytical 
solution of the Stokes equation with boundary conditions Eq. f[T2j) at the wavy wall with am- 
plitude a and wavelength A was obtained for two limiting cases and small ka. For Lq/X ^ 1 
the effective slip length is given by 

lim L e jj= Lq — ka 2 u (ka), (14) 

while in the limit of Lq/\ ^> 1 it reduces to 

v ( 1 k3(l2 V 1 

hm L eff = — + — - , 15 

feLo^oo VL uj^ika)/ 

where the functions u (ka) and uj^ka) are defined as 

_ 1 - l/4(fea) 2 + 19/64(fcq) 4 + Q[(kaf] 
^ ka >- l + (k)2_ 1/2N 4 + 0[(k) 6] > ( lb > 

and 

_ 1 - 5/4(fca) 2 + 61/64(fcq) 4 + 0[{kaf] 
UJ ^ ka >~ l + {kaf-l/2{kaY + 0[{kaf} ' { } 

An approximate analytical expression for the effective slip length that interpolates between 
the two bounds Eq. ( fl4l) and Eq. (115]) is given by 

Louj^ka) - ka 2 LU (ka)/(l + 2kL ) 
Leff ~ 1 + k 3 a 2 L ' (18) 
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with the range of applicability ka < 0.5. For larger wavenumbers ka > 0.5, the function 
Uoo^ka) overestimates the numerical solution and the interpolated formula Eq. ffl8|) does not 



B. Flow over a rough surface with the local no-slip boundary condition 

In this section, the Stokes equation with the local no-slip condition at the upper and lower 
walls is solved numerically to study the effect of corrugation amplitude on the effective slip 
length. The velocity profiles, averaged over the period of corrugation A, are plotted in the 
inset of Fig.[3]for several values of wavenumber ka = 2ira/ A. As ka increases, the normal- 
ized velocity profiles remain linear in the bulk region and become curved near the lower 
corrugated wall. The linear part of the velocity profiles is used to compute the effective slip 
length, which is plotted as a function of wavenumber in Fig.[3j With increasing corruga- 
tion amplitude of the lower wall, the effective slip length decays monotonically and becomes 
negative, indicating that the effective no-slip boundary is shifted into the fluid domain. For 
ka < 1, the numerical results agree well with the analytical solution Eq. (fl8|) denoted by 
the solid line in Fig.[3j The deviation from the analytical solution becomes significant at 
larger wavenumbers where the streamlines extracted from the Stokes solution indicate the 
presence of backflow at the bottom of the valley. 

In order to investigate the flow behavior above the sinusoidal surface, the shear stress 
and pressure along the lower wall were computed from the solution of the Stokes equation. 
In the presence of surface curvature the wall shear stress r w has two components 



where dut/dn is the normal derivative of the tangential velocity Ut, and R(x) is the local 
radius of curvature. In the case of no-slip boundary condition {u t = 0), the local shear stress 
at the wall is reduced to r w = The normalized shear stress along the lower corrugated 

wall is plotted in Fig.H] for different corrugation amplitudes. The maximum value of the 
shear stress is located at the peak of the surface corrugation (x/X = 0.25) and it increases 
with increasing amplitude, which is consistent with the results of previous analytical studies 
of a laminar flow over a wavy wall 0, |2| . The fluid tangential velocity near the boundary 
is proportional to the wall shear stress shown in Fig.Hl Therefore, the tangential velocity 



apply [26 1 . 




(19) 



S 



is also maximum above the peak and, as the flow moves downstream, it decelerates and 
the velocity becomes zero inside the valley at sufficiently large amplitudes. For ka ^ 0.79, 
the shear stress profiles intersect the dashed line (r w = 0) at two points and a clockwise 
flow circulation develops inside the valley. As the corrugation amplitude increases, the 
intersection points move away from each other and the flow recirculation region becomes 
larger. These results are in agreement with previous estimates of the critical wavenumber 
ka « 0.77 for the onset of flow separation in sufficiently thick films 44], 



The pressure along the lower wavy wall is plotted in Fig.[5]for the same amplitudes as 
in Fig.Hl The value P* used for normalization is the maximum pressure, which is located 
above the wavy surface with ka = 1.12 on the left side of the peak. For each wavenumber, 
the pressure along the surface is maximum on the left side of the peak, where the surface 
faces the mainstream flow. The surface pressure reaches its minimum value on the right side 
of the peak (see Fig.|3j). As the flow moves further downstream into the valley, it encounters 
an adverse pressure gradient, which becomes larger as ka increases. At large wavenumbers 
ka ^ 0.83, the flow near the surface cannot overcome the combined resistance of the viscous 
forces and the adverse pressure gradient, and separates from the surface at the point where 

T w = 0. 

The pressure contours and streamlines near the corrugated surface with wavenumber 
ka = 1.12 are depicted in Fig. 0(a). The pressure contours indicate the presence of an 
adverse pressure gradient in the region 0.3 < x/X < 0.6 on the right side of the peak (see 
also Fig.EJ). The streamlines illustrate the flow separation inside the valley at x/X ~ 0.52. 
After the separation point the flow near the wall reverses direction and moves against the 
mainstream. The local velocity profile inside the valley is shown in the inset of Fig.[6](a). 
As the flow cross-section decreases towards the peak, the flow attaches to the surface at 
the point where the shear stress is zero again. Note also that for ka ^ 0.83 in Fig. [5j the 
pressure profiles along the lower surface exhibit a slight drop inside the valley where the 
vortex is present. After the attachment point, due to the periodic boundary condition in 
the x direction, the surface pressure increases up to its maximum value, which is located on 
the left side of the peak. The nonlinearity due to inertia is absent in the Stokes flow and 
the vortex in the valley remains symmetrical. 

We also comment that a secondary vortex is formed inside a deeper cavity (ka « 2.33), 
which is counter-rotating with respect to the primary vortex (not shown). Previous an- 
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alytical study of a creeping flow over a wavy wall demonstrated that a secondary vortex 
appears at ka > 2.28 when the channel width is larger than the corrugation wavelength 45]. 
In the limiting case, when the cavity consists of two parallel walls, an infinite sequence of 



counter-rotating vortices appears between the walls 
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C. Effect of the local slip on the flow pattern near the rough surface 

Next, we present the results of the numerical solution of the Stokes equation with the 
local slip condition at the lower wavy wall while the boundary condition at the upper wall 
remains no-slip. The pressure contours and the streamlines are plotted in Fig.[6]for several 
values of the intrinsic slip length Lq and the wavenumber ka = 1.12. The streamline patterns 
indicate that with increasing slip length Lq, the size of vortex inside the valley is progressively 
reduced and the vortex eventually disappears. Similar to the analysis in the previous section, 
the pressure and shear stress along the lower wall are computed in the presence of the local 
boundary slip. 

The normalized pressure along the lower wavy wall is plotted in Fig. [7] as a function of 
the slip length Lq. Similar to the case of no-slip boundary condition, the profiles exhibit a 
maximum and a minimum in pressure on the right and left sides on the peak, respectively. In 
a wide range of L , the difference between the locations and the magnitudes of the extrema 
is barely noticeable. As the flow moves down a slope (0.3 < x/X < 0.5 region in Fig.[7j), 
the pressure gradient along the surface becomes positive and its magnitude decreases with 
increasing values of Lq. The separation and attachment points are located at the intersection 
of the shear stress profiles with the horizontal line (t w = 0) shown in Fig.[HJ In comparison 
to the no-slip case, the smaller combined effect of the adverse pressure gradient and the 
wall shear stress causes a shift of the separation point deeper into the valley (e.g. see 
Fig. EI). As the slip length Lq increases, the separation and attachment points move closer 
to each other, the vortex becomes smaller and eventually vanishes. We also note that at the 
separation point both components of the shear stress du t /dn and u t /R{x) become zero in 
agreement with the MRS criterion for the flow separation of the boundary layer at a moving 



substrate [52] 



The influence of the local slip on the shape of the velocity profile and the normal derivative 
of the tangential velocity at the bottom of the valley is illustrated in the insets of Fig. El 
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The local slip length L s (x) in Eq. (|T2|) is a function of the radius of curvature R(x) and 
the intrinsic slip length Lq. For all cases considered in Fig. [6], the condition Lq < \R{x)\ 
holds and the local slip length L s remains positive everywhere along the lower wall, which 
means that the slip velocity u t and the normal derivative du t /dn carry the same sign. It is 
expected, however, that for the opposite condition, L > |-R(x)|, the values u t and du t /dn 
would have different signs at the bottom of the valley, and the corresponding velocity profile 
will be qualitatively similar to the profile shown in the inset of Fig.[6]^d) but shifted by a 
negative slip velocity (see next section). 

The effective slip length computed from the numerical solution of the Stokes equation and 
Eq. (I18p is plotted in Fig. as a function of the intrinsic slip length Lq for wavenumbers ka = 
0.28 (a/h = 0.04) and ka = 1.12 (a/h = 0.18). For small values of L Q , the effective slip length 
approaches a negative value previously reported in Fig.[3]for the no-slip case. As Lo increases, 
L e jj grows monotonically and appears to saturate to a constant value. The transition of the 
effective slip length from a growing function of Lq to a nearly constant value occurs at 
larger Lq when ka decreases. For the small wavenumber ka = 0.28, the effective slip length 
computed from Eq. ( {TBI) is in a good agreement with the numerical solution of the Stokes 
equation. Visual inspection of the streamline patterns indicates that there is no backflow at 
any Lq. In the saturation regime, Lq — > oo, the wall shear stress becomes zero everywhere 
along the lower boundary, the streamlines near the lower wall follow the boundary curvature, 
and the effective slip length in Eq. ( fT5l) approaches L e ff/X ~ 1/ [27i(ka) 2 ] — 9/8ir. For the 
large wavenumber ka = 1.12, the analytical results Eq. ( |T8l) overestimate L e g computed from 
the numerical solution at Lq/Ii > 0.02 and the flow circulation is developed in the valley at 
L /h ^ 0.067. The vortex vanishes at the bottom of the valley at sufficiently large values of 
Lq (denoted by the vertical arrow in Fig.[9]), and the flow streamlines are deformed to follow 
the boundary curvature [e.g. see Fig.[6]^d)]. The results for the intrinsic slip length, which 
determines the threshold for the onset of the flow circulation at the bottom of the groove, 
are summarized in Fig.[TQJ For the wavenumbers examined in this study, ka ^ 1.26, the 
numerical simulations indicate that if the flow circulation is present in the valley then the 
effective slip length is negative and L e g increases with decreasing vortex size. 
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D. Effect of the Reynolds number on the effective slip length 

The analysis of the Stokes equation discussed in the previous section demonstrated that 
increasingly large local slip at the lower wavy wall eliminates the flow circulation in the 
valley and leads to a larger effective slip. In this section, the influence of the inertia term in 
the Navier-Stokes equation on the flow pattern and the effective slip length is investigated. 
For the shear flow with slip condition at the lower corrugated wall L , the Reynolds number 
is defined as 

Re = pg-m+Wft), (20) 

where p is the fluid density, U* is the upper wall velocity, and h (1 + L e jj/h) is the distance 
between the upper flat wall and the effective no-slip boundary plane. In the case of no-slip 
boundary condition at the lower flat wall, L e jj is zero and the standard definition of the 
Reynolds number is recovered, i.e., Re = pU*h/ p. 

The pressure and shear stress along the lower wavy wall are plotted in Fig. [11] for the 
selected values of the Reynolds number and no-slip boundary conditions. As Re increases, 
the adverse pressure gradient along the right side of the corrugation peak (0.3 < x/X < 0.5) 
becomes larger and the pressure drop inside the valley (0.5 < x/X < 0.9) increases. For 
each value of the Reynolds number in Fig. fTTT b). similar to the Stokes flow case, the shear 
stress above the peak is maximum and, as the flow moves downstream along the right side 
of the peak, it decelerates and eventually separates from the surface when r w = 0. With 
increasing upper wall velocity, the shear stress above the peak increases and causes the 
flow to decelerate faster along the right side of the peak. The separation and attachment 
points (determined from the condition r w = 0) move further apart from each other and the 
circulation region inside the valley expands [see Fig. fTlT b)]. 

The effect of the inertia term in the Navier-Stokes equation on the shape of the wall shear 
stress and pressure profiles can be seen in Fig.[TTJ The shear stress profiles above the peak 
are not symmetric with respect to x/X = 0.25, which means that the flow decelerates faster 
on the right side of the corrugation peak than it accelerates on the left side. Also, the average 
adverse pressure gradient on the left side of the peak is larger than its value on the right 
side of the peak (0.3 < x/X < 0.5). By increasing the upper wall velocity, the separation 
point moves further upstream than the attachment point downstream. The formation of 
asymmetric vortex flow at finite Reynolds numbers is consistent with previous findings for a 
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flow in an undulated tube |36|, |39_|, |40|, |42( . The pressure contours and streamlines extracted 
from the NS equation with the no-slip boundary condition are plotted in Fig. [T27 a) for 
ka = 1.12 and Re = 79. The flow streamlines in the valley indicate an asymmetric clockwise 
circulation, which is larger than the flow circulation region shown in Fig. 0(a) for the Stokes 
case. 

In the presence of the local slip condition along the lower corrugated wall, the size of 
the vortex becomes smaller while the flow structure remains asymmetric (see Fig. [T2|) . The 
decrease of the vortex size is similar to the case of the Stokes flow shown in Fig. [6] and can 
also be described in terms of the pressure and shear stress along the lower wall. In Fig.fTBl 
the pressure and shear stress profiles are plotted for the same values of the upper wall 
velocity as in Fig. [TT] but with the slip boundary condition (L /h = 0.25) along the lower 
wall. Note that the Reynolds numbers in Fig. [13] are slightly larger than the values reported 
in Fig. [11] for the same U* because of the larger effective slip length entering the definition 
of the Reynolds number [see Eq. ( [20]) ]. For each value of the upper wall velocity, the adverse 
pressure gradient and the wall shear stress on the right side of the peak are smaller than 
in the case of the no-slip boundary condition, and, as a result, the vortex either becomes 
smaller or vanishes (see Figs.fTT] and [TSl . As the Reynolds number increases, however, the 
vortex forms and expands asymmetrically to fill the bottom of the groove. The inset of the 
Fig.fTWb) demonstrates that the slip velocity at the bottom of the valley is negative while 
its normal derivative dut/dn is positive, in contrast to the velocity profiles shown in Fig.[6] 
for the same ka = 1.12 and smaller slip lengths Lo/h ^ 0.08. 

The effective slip length is plotted in Fig. [14] as a function of the Reynolds number for 
the selected values of the intrinsic slip length L and ka = 1.12. With increasing Reynolds 
number, the flow streamlines move away from the lower boundary and straighten out, the 
slope of the normalized velocity profiles increases and the effective no-slip boundary plane 
is shifted into the fluid domain. For Lo/h ^ 0.067, the circulation is always present in the 
valley and the flow streamlines in the bulk of the fluid do not penetrate deep into the valley 
[e.g. see Fig. [T27 b)]. For larger slip lengths Lo/h > 0.067, the flow streamlines show that 
there is no backflow at low Re, and the vortex is formed at the bottom of the valley only at 
sufficiently large Reynolds numbers indicated by the dashed line in Fig.[T31 The numerical 
results obtained from the solution of the Navier-Stokes equation demonstrate that the growth 
or decay of the vortex as a function of the Reynolds number or the intrinsic slip length is 
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accompanied by the decrease or increase of the effective slip length. 

IV. CONCLUSIONS 

In this paper the effects of local slip boundary condition and the Reynolds number on 
the flow structure near sinusoidally corrugated surfaces and the effective slip length were 
investigated numerically by solving the Stokes and Navier-Stokes equations. The effective 
slip length was defined with respect to the mean height of the surface roughness by extrap- 
olating the linear part of the velocity profile averaged over the corrugation period. In the 
case of the Stokes flow with the local no-slip boundary condition, the effective slip length 
decreases with increasing corrugation amplitude and the vortex flow develops in the groove 
of the rough surface for ka ^ 0.79. In the presence of the local slip boundary condition 
along the wavy wall, the effective slip length increases and the size of the recirculation zone 
is reduced. The vortex vanishes at sufficiently large values of the intrinsic slip length. The 
analysis of the pressure and wall shear stress computed from the Navier-Stokes equation 
shows that the asymmetric vortex flow develops in the groove due to the inertia term even 
when the local slip boundary condition is applied. The effective slip length decreases with 
increasing Reynolds number. The numerical simulations suggest that the variation of the 
vortex size as a function of either the Reynolds number or the intrinsic slip length correlates 
with the magnitude of the effective slip length. 
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FIG. 1: (Color online) Schematic diagram of the steady-state Couette flow over a rough surface. 
The upper flat wall is moving with a constant velocity U in the x direction. The lower station- 
ary wall is modeled as a sinusoidal wave with amplitude a and wavelength A. The wavenumber 
ka = 2na/X varies in the range ^ ka ^ 1.26. 
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FIG. 2: Diagram of a bilinear element in (a) the physical coordinate system (x, z) and (b) a 
transformed element in the natural coordinate system (£,77). 
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FIG. 3: (Color online) The effective slip length as a function of wavenumber ka computed from 
the solution of the Stokes equation with the no-slip boundary condition. The primary vortex is 
formed at the bottom of the valley for ka ^ 0.79 (see the vertical arrow at ka = 0.79). The solid 
line is calculated using Eq. (|18p . The inset shows the normalized velocity profiles obtained from 
the Stokes solution for the selected values of ka. The dashed line located at a = is the reference 
for calculating the effective slip length. 
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FIG. 4: (Color online) Shear stress along the lower wavy wall computed from the Stokes solution 
with no-slip boundary condition (Lq = 0) for the indicated values of wavenumber ka. The value 
r* used for normalization is the shear stress at the flat wall. The intersection of the curves with 
the dashed line determines the location of the flow separation and attachment inside the valley. 
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FIG. 5: (Color online) The normalized pressure along the lower corrugated boundary extracted 
from the Stokes solution with no-slip boundary condition for the same values of ka as in Fig. [31 
The maximum pressure P* for ka = 1.12 and Lq = is located in the bulk region at (x/X, z/X) ~ 
(0.07,0.23). 
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FIG. 6: (Color online) Pressure contours and streamlines for the wavenumber ka = 1.12 and the 
slip length at the lower wall L = (a), L /h = 0.03 (b), L /h = 0.06 (c), and L /h = 0.08 (d). 
The pressure contours are normalized by the maximum value P* located on the left side of the 
peak (x/X,z/X) ~ (0.07,0.23) in the case of Lq = 0. The vertical dashed line inside the valley at 
x/X = 0.75 indicates the cross-section used to compute the velocity profiles shown in the inset. 
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FIG. 7: (Color online) Normalized pressure along the lower wavy wall extracted from the solution 
of the Stokes equation for the tabulated values of Lq and ka = 1.12. The maximum pressure P* 
for Lq = is located above the lower surface at (x/X,z/X) ~ (0.07,0.23). 
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FIG. 8: (Color online) Shear stress along the corrugated lower wall (ka = 1.12) for the indicated 
values of the local slip length Lq. The intersection of the dashed line with the shear stress profiles 
shows the location of the flow separation and attachment inside the valley. 
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FIG. 9: The effective slip length, extracted from the Stokes solution as a function of Lq for ka = 0.28 
(top) and ka = 1.12 (bottom). The dashed line is computed from Eq. (|18p . For ka = 1.12 and 
Lo/h ^ 0.067, the vortex is formed in the valley (see the vertical arrow at L^/h = 0.067). The 
error bars are smaller than the symbol size. 
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FIG. 10: (Color online) The intrinsic slip length Lq above which there is no vortex at the bottom 
of the valley (top) and the corresponding effective slip length L e g (bottom) computed from the 
Stokes solution and plotted as a function of wavenumber ka. The increment in the slip length used 
to determine the threshold values is about the symbol size (top). The dashed lines are quadratic 
fits to guide the eye. 
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FIG. 11: (Color online) The normalized pressure (top) and shear stress (bottom) along the lower 
wavy wall as a function of the Reynolds number for ka = 1.12 and no-slip boundary conditions. 
The value P* is the maximum pressure located above the lower boundary on the left side of the 
peak (x/X, z/X) ~ (0.07, 0.23) for L = and Re = 4.0. 
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FIG. 12: (Color online) Pressure contours and streamlines extracted from the solution of the 
Navier-Stokes equation for ka = 1.12 and Lq = (a) and Lq/H = 0.25 (b). The value P* is the 
maximum pressure located at (x/A, z/X) ~ (0.07,0.23) for ka = 1.12, Lq = 0, and Re = 4.0. The 
dashed line inside the valley indicates the cross-section used to compute the velocity profiles shown 
in the inset. 
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FIG. 13: (Color online) The normalized pressure (top) and shear stress (bottom) along the lower 
wavy wall as a function of the Reynolds number for ka = 1.12 and Lq/H = 0.25. The normalization 
value P* is the maximum pressure located above the surface on the left side of the peak (x/X, z/X) ~ 
(0.07,0.23) for L = and Re = 4.0. 
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FIG. 14: (Color online) The effective slip length computed from the solution of the Navier-Stokes 
equation as a function of the Reynolds number for ka = 1.12 and Lq = (>), Lo/h = 0.09 (a), 
L0//1 = 0.12 (<), and Lo//t = 0.25 (▼). The dashed line indicates the upper bound of the region 
where a vortex is present in the groove of a rough surface. The error bars associated with the 
threshold values of the Reynolds number are about the symbol size. 
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